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ABSTRACT 



Context. Deep photometric surveys for substellar objects in young clusters and for high-redshift quasars are affected by contaminant 
sources at different heliocentric distances. If not correctly taken into account, the contamination may have a strong effect on the Initial 
Mass Function determination and on the identification of quasars. 

Aims. We calculate in detail the back- and foreground contamination by field dwarfs of very late spectral types (intermediate and late 
M, L and T) in deep surveys and provide the data and tools for the computation. 

Methods. Up-to-date models and data from the literature have been used: (/) a model of the Galactic thin disc by an exponential law; 
((7) the length and height scales for late-type dwarfs; (Hi) the local spatial densities, absolute magnitudes and colours of dwarfs for 
each spectral type. 

Results. We derive a simplified expression for the spatial density in the thin disc that depends on the heliocentric distance and the 
galactic coordinates (/, b) and integrate it in the truncated cone screened in the survey. As a practical application, we compute the 
numbers of L- and T-type field dwarfs in very deep (/ = 21-29 mag) surveys in the direction of the young a Orionis cluster. The 
increasing number of contaminants at the faintest magnitudes could inhibit the study of the opacity mass limit at M < 0.003 M Q in 
the cluster. 

Key words, stars: low mass, brown dwarfs - stars: luminosity function, mass function - Galaxy: stellar content - open clusters and 
associations: individual: a Orionis - methods: analytical 



1. Introduction 

Deep wide-field multi-band photometric searches in open clus- 
ters and star-forming regions are the most useful tool for the dis- 
covery and identification of young brown dwarfs and planetary- 
mass objects (e.g., Rebolo et al. 1996 and Bouvier et al. 1998 
in the Pleiades; Preibisch & Zinnecker 1999 and Ardila, Martin 
& Basri 2000 in Upper Scorpius; Zapatero Osorio et al. 2000 
and Bejar et al. 2001 in cr Orionis). Unfortunately, an impor- 
tant fraction of the selected cluster member candidates do not 
belong to the cluster, but actually are interlopers in the fore- 
and background. This fact strongly affects posterior analyses of 
the bottom of the Initial Mass Function, down to about a few 
Jupiter masses. There are several ways of reducing the contam- 
ination of the mass function incurred by these interlopers, in- 
cluding the identification of spectral indicators of youth/low sur- 
face gravity in individual sources (Martin, Rebolo & Zapatero 
Osorio 1996; Stauffer, Schultz & Kirkpatrick 1998; White & 
Basri 2003; McGovern et al. 2004; Mohanty, Jayawardhana & 
Basri 2005) and removing sources with kinematics inconsistent 
with the cluster (Stauffer, Hamilton & Probst 1994; Hambly et al. 
1999; Moraux, Bouvier & Stauffer 2001; Bihain et al. 2006; 
Lodieu et al. 2007b). However, these techniques may require 
prohibitive investments of time on large telescopes, particularly 
for the lowest mass objects (e.g., S Ori 70 - see below), and 
may not eliminate contaminants for distant and/or older clus- 
ters and associations. For example, at 100 pc, accuracies better 
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than ~10mas/a _I are needed to measure velocities of ~5kms~ 1 . 
In that case, there can still be foreground objects that share the 
same motions but are not cluster members. Their membership 
status can only be insured through the detection of youth signa- 
tures. Some known features of extreme youth (age < 10 Ma) are 
Li i /16707.8 A in absorption (Rebolo, Martin & Magazzu 1992), 
abnormal strength of alkali lines and water and H2 absorption 
bands due to low gravity (Martin et al. 1999; Kirkpatrick et al. 
2006; Allers et al. 2007 - H2 can also be in emission from 
discs [Thi et al. 2001]), strong Ha /16562.8A emission asso- 
ciated to accretion and mid-infrared flux excesses due to cir- 
cum(sub)stellar discs. However, obtaining high signal-to-noise 
ratio, mid-resolution spectroscopy or Spitzer Space Telescope 
imaging of the faintest, reddest objects in a deep photometric 
survey (/ > 18-19mag, J > 15-16mag) can be a very diffi- 
cult task. 

A useful statistical approach is to measure both the luminos- 
ity function in the area of the cluster and in neighbouring fields, 
and extracting the true cluster luminosity function as the dif- 
ference of these. This method has been succesfully employed 
in the Orion Nebula Cluster and described in detail by, e.g., 
Muench et al. (2002). In addition to relying on a local extinc- 
tion map of the cluster to redden the stellar population in the 
extinction-free comparison fields, the method has, regrettably, 
an important handicap: one must integrate the same exposure 
time in the neighbouring fields to properly subtract the "field 
mass function". For the deepest surveys in clusters, that use ex- 
pensive large facilities, this time-consuming solution cannot be 
accomplished in many cases. As a result, the vast majority of 
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these deep surveys establish certain selection criteria for cluster 
membership based on colour-magnitude diagrams. 

Depending on the criteria and the number and suitability 
of passbands, more or less cluster non-member contaminants 
could be among the selected member candidates in a photo- 
metric survey. To define a selection criterion, the authors must 
account for photometric uncertainties in the survey, the natural 
scatter of the cluster sequence (i.e., the sources themselves do 
not have a single absolute magnitude or colour, but rather a range 
of both appropriate for their spectral type), unresolved binarity 
and intrinsic photometric variability. Figure 2 in Caballero et al. 
(2007) illustrates a typical selection of member candidates in 
the cr Orionis cluster based on the loci of spectroscopically con- 
firmed members in an / vs. /- J diagram. Theoretical isochrones 
may be used as reference, as well. The removal of interlop- 
ers, whose location coincides with that of cluster members in 
colour-magnitude diagrams, is very important in the faintest and 
reddest magnitude intervals. Quantitative estimations of contam- 
inants have been presented in several searches in young clus- 
ters (Paresce, de Marchi & Romaniello 1995; Zapatero Osorio, 
Rebolo & Martin 1997; Martin et al. 2000; Lucas et al. 2001; 
Jeffries et al. 2004; Moraux et al. 2007). 

In this paper we provide the necessary tools and state-of-the- 
art data for the correct decontamination of field late-type stars 
and brown dwarfs in deep photometric surveys at intermediate 
and high galactic latitudes. A practical application of the inte- 
grated number of contaminants in very deep surveys towards 
the o~ Orionis cluster is presented. Our work can be applied to 
searches in other open clusters and star-forming regions, but also 
to high-z quasar surveys. 

2. Analysis 

2. 1 . Possible contaminants 

The contaminants appearing in deep photometric searches for 
very red objects are in general faint reddened galaxies, variable 
subgiants, (reddened) distant giants, or field dwarfs with very 
late spectral types that fall close to the sequence of the cluster 
under study in a colour-magnitude diagram: 

- Galaxies are in general extended sources, so their FWHMs 
are larger than the average FWHM of the stars in the studied 
field of view, and they are easily rejected during the point- 
spread-function photometry: down to J ~ 19-20 mag, galax- 
ies can be revealed under good seeing conditions (see, for 
example, the recent work by Foster et al. [2008]). They also 
display colours that do not match the dwarf sequence in a 
colour-colour diagram (e.g., J - K s vs. / - J). Some re- 
mote quasars with z ~ 6 have, however, point-like FWHMs 
and spectral energy distributions nearly identical of those of 
mid-T dwarfs; these quasars are, nevertheless, extremely rare 
(Stern et al. 2007). 

- Except for extremely young star forming regions with high 
extinction clouds (p Ophiuchi, Chamaeleon I, Orion Nebula 
Cluster), open clusters are relatively free of interstellar ab- 
sorption. In some clusters, as cr Orionis, the strong wind and 
ultraviolet radiation from the Trapezium-like central system 
has completely cleared the region from cometary globulae, 
molecular clouds and most of the original intracluster ma- 
terial. In the case of clusters with thick molecular clouds, 
these can act as a "shield to background field stars" (Lucas 
& Roche 2000; Luhman et al. 2000; Muench et al. 2002), 
which may facilitate the decontamination (there will be only 
interloping non-cluster field stars in foreground). 



- According to Kirkpatrick et al. (1994), contamination by gi- 
ant stars in areas excluding the Galactic plane is negligi- 
ble. Given the absolute magnitude of late-type giant stars, 
a hypothetical faint (intrinsically bright) giant contaminant 
would be at several kiloparsecs over or below the Galactic 
plane, where only the globular clusters are found. More re- 
cent works have shown a comparatively large number of gi- 
ant interlopers. For example, Cruz et al. (2007) identified 
33 giants and 17 carbon stars among a list of ~200 ultra- 
cool dwarf candidates in the field. However, they only used 
2MASS JHK S photometry for the selection. The employ- 
ment of at least one optical filter and astrometric and spec- 
troscopic information from the literature can help eliminate 
these contaminants in cluster (and field) surveys. 

- Photometric variability of red (sub-)giants in the background 
is an additional source of contamination for surveys whose 
multiband images were taken at different epochs (i.e., the 
colours of variable sources are not "real"). Many of such 
contaminants, although difficult to discard a apriori, can be 
identified by studying their light curves. Very often, the deep 
final image in a pass band of a survey is the combination of 
a long series of exposures (this methodology, if accompa- 
nied with a suitable dithering, allows to reach deeper magni- 
tudes than in a single exposure of the same total time). See 
Caballero et al. (2004) and Scholz & Eisloffel (2004) for ex- 
amples of photometric variability studies in a young open 
cluster based on such long series. 

- The major contributors to red contamination in deep surveys 
in clusters are, therefore, field very low-mass stars and brown 
dwarfs with late spectral types. The discovery of distant, late- 
type dwarfs in deep surveys is not exceptional. For exam- 
ple, Stanway et al. (2008) found a population of M dwarfs 
at large heliocentric distances in the GOODS fields. Pirzkal 
et al. (2005) also identified 18 M- and two L-type candidate 
dwarfs in the Hubble Ultra Deep Field, that require further 
observations for confirmation. Many more field L dwarfs 
have been found in the foregrounds of young open clus- 
ters and associations (e.g., cr Orionis: Caballero et al. 2007; 
Upper Scorpius: Lodieu et al. 2008). Cuby et al. (1999) iden- 
tified during an extragalactic survey the most distant field 
T dwarf detected to date (NTTDF J 1205-0744), at an he- 
liocentric distance d ~ 90 pc, and with a magnitude J ~ 
20 mag. Kendall et al. (2007) and Lodieu et al. (2007a) re- 
cently estimated slightly lower distances (d ~ 50-80) for 
the intermediate T dwarfs ULAS J145243.59+065542.9 and 
ULAS J223955.76+003252.6, respectively. The actual sta- 
tus of SOri 70, a ~T6-type, ~3Mj up , planetary-mass ob- 
ject candidate towards Orion, is still under debate (Zapatero 
Osorio et al. 2002b, 2008; Martm & Zapatero Osorio 2003; 
Burgasser et al. 2004; Scholz & Jayawardhana 2008). 



2.2. The Galaxy thin-disc model 

Within the context of the "standard" Galactic model, the Galaxy 
can be modelled by a double exponential, that represents the 
thin and the thick discs, and a power law for the halo (e.g., a 
de Vacouleurs-type power law). Classic works and reviews on 
Galactic structure were presented by, e.g., Bahcall & Soneira 
(1980), Gilmore & Reid (1983), Majewski (1993) and Kroupa, 
Tout & Gilmore (1993). Recent observational works on Galaxy 
disc models from deep imaging surveys can be found in, for ex- 
ample, Phleps et al. (2005), Ryan et al. (2005), Karaali (2006), 
and Juric et al. (2008). 
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The volume densities of the axisymmetric thin and thick 
discs depend on the galactocentric distance R and on the height 
over or below the galactic plane, Z. For an object of spectral type 
i, the expression for its spatial density depending on heliocentric 
galactic coordinates (I, b) and heliocentric distance (d) is usually 
written as follows: 

R{d,lM-Ro [Z +d sin b\ 

rij(d,l,b) = noje e h z (1) 

(e.g., Chen et al. [2001] and references therein), where no,- is the 
star spatial density at the Galactic plane (Z = 0) and at the solar 
galactocentric distance (R = R G ), Z is the height of the Sun 
over the Galactic plane, and Iir and hz are the length (radial) and 
height scales, respectively. The exponential vertical distribution 
actually comes from the approximation sech 2 (|Z|/2/iz) ~ e^ z ^ hz 
(Gilmore, King & van der Kruit 1990). We discard using more 
complex forms for the vertical density profile (e.g., Zheng et al. 
2004), since additional parameters that characterize these forms 
(e.g., yS) are more poorly constrained than the exponential scale 
height factors (see also the following discussion on the thin-to- 
thick disc normalisation). The galactocentric distance R(d, I, b) 
of the object is, in its turn: 

R(d, I, b) = {R% + d 2 cos 2 b - 2R dcos b cos l) 1/2 . (2) 

(e.g., Bahcall & Soneira 1980 - R is measured along the Galactic 
plane). Assuming that the heliocentric distance, d, is much less 
than the solar galactocentric distance, d <k R q (see below for a 
justification), then the galactocentric distance of the object can 
be approximated by: 

R(d,l,b) *R Q -dcosbcosl. (3) 

The minus sign guarantees that R < R G when the line of sight 
points to the Galactic Centre (I = deg), and vice versa (we are 
using the right-handed Cartesian galactocentic coordinate sys- 
tem in which the X axis points towards the Sun, the Y axis in 
the anti-direction of Galactic rotation and the Z axis towards the 
North Galactic Pole). In the linear approximation, Eq.[TJcan be 
rewritten as: 

rii(d,l,b) « noje e ;, z = n^e e '•« h z , (4) 

where the sign convention depends on whether d sin b is greater 
or less than -Z Q (i.e., whether the source is above or below the 
Galactic plane). Defining two auxiliar variables that vary with 
the galactic coordinates, riA,i an d dgih b), Eq.[T]results on a sim- 
ple expression with the heliocentric distance: 

rn(d,l,b) * riAje' 1 ^, (5) 
where the auxiliar variables are computed through: 

riAj = n je h ? (6) 
and 

1 cos b cos/ sinb 
= + . (7) 

d B (l, b) h R h z 

In Table Q] we show the main parameters of the Galaxy thin- 
disc model from Chen et al. (2001). These values suitably fit 
classic and recent determinations (e.g., Faber et al. 1976; Bahcall 
& Soneira 1984; Kerr & Lynden-Bell 1986; Kuijken & Gilmore 
1989; Robin, Reyle & Creze 2000; Pirzkal et al. 2005; Juric et al. 



Table 1. Main parameters of the Galaxy thin-disc model from 
Chen et al. (2001). 



Parameter 


Re 




h R 


hz 




(pc) 


(pc) 


(pc) 


(pc) 


Value 


8 600+200 


+27±4 


2 250+1000 


330±3 



2008 and references above). We will subsequently use only Chen 
et al.'s values for consistency, even though there are better es- 
timates for some parameters (e.g., the solar galactocentric dis- 
tance - Reid 1993; Eisenhauer et al. 2003). From the constancy 
of the height scale of GKM-type dwarfs, as shown by Bahcall & 
Soneira (1980) and Gilmore & Reid (1983), it seems reasonable 
to also use the same height scale for M, L and T dwarfs. Pirzkal 
et al. (2005) measured h z = 400+100 pc for late M and L dwarfs, 
which is a little larger than ours, but still consistent. While the 
vertical scale heights of early-type stars (O, B, A) are shorter, 
they will not typically contribute to the photometrically selected 
sample and are rare anyways. 

The thick disc is more rarified and extended than the thin 
disc, with a thin-to-thick disc normalisation factor in the so- 
lar vicinity of 13-6.5% (Chen et al. 2001; Juric et al. 2008). 
This value interval is larger than previously reported, like the 
2% in Gilmore (1984) or the lack evidence for a thick disc 
in Bahcall & Soneira (1984). According to theoretical cooling 
sequences (Burrows et al. 1997; Chabrier et al. 2000; Baraffe 
et al. 2003), very old thick-disc substellar objects with very 
low masses (M <sc 0.08 M Q ) may be dimmed to exceedingly 
faint magnitudes, which would prevent their detection (thick- 
disc stars and brown dwarfs are thought to be very old, with 
ages > 10 Ga - e.g., Fuhrmann 1998; Prochaska et al. 2000; 
Feltzing, Bensby & Lundstrom 2003). For example, an 0.05 Mo- 
mass brown dwarf dims from M/ = 17.2 mag at 1 Ga to M/ = 
20.5 mag at 10 Ga. The same brown dwarf at that time intervals 
would have effective temperatures typical of L5-7V and T7- 
9V, respectively. Since very low-mass stars above the hydrogen 
burning limit (with late M and very early L spectral types) keep 
a roughly constant luminosity in the main sequence for tens gi- 
gayears, we expect that the ratio between L and T (and cooler) 
dwarfs in the (relatively old) thick disc is less than in the (rela- 
tively young) thin disc. In other words, field late-M-, L- and T- 
type dwarfs of the solar neighbourhood are younger, in average, 
than G-, K- and early M-type dwarfs (see Table 5 in Zapatero 
Osorio et al. 2007a). There is a discussion in the negligibility of 
thick disk and halo brown dwarfs in a given imaging sample in 
Pirzkal et al. (2005). Both the low thin-to-thick disc normalisa- 
tion factor and the relative faintness of thick-disc objects suggest 
that it is appropriate to use just the thin-disc exponential for de- 
termining the density of late-type dwarfs in the Galaxy. The error 
in this assumption (< 10 %) is smaller than, or of the order of, the 
uncertainties in the determination of R G , Z and Iir. 

Using the values in Table[T| the ratio riA,i/no,i = e *z (Eq. |6J) 
equals 0.921+0.011 for Z > 0, and 1.085+0.013 for Z < and 
is independent of the galactic coordinates. However, the auxiliar 
variable dsiU b) strongly depends on the galactic longitude and 
latitude in the direction of observation and has a large interval 
of variation. 

Fig. [T] shows the suitability and accuracy of the linear ap- 
proximation of R(d, I, b). The differences between the relative 



4 



J. A. Caballero et al.: Field late-M, L and T dwarfs in deep surveys 




Table 2. Heliocentric distances and galactic coordinates of three 
representative young open clusters. 



d[pc] 

Fig. 1. Top window: variation of the relative spatial density of 
stars and brown dwarfs of spectral type i in the thin disc with 
the heliocentric distance from 5 to lOOOOpc in the direction 
of three representative young clusters: Upper Scorpius (dotted 
line), Pleiades (dash-dotted line), and <x Orionis (-red- solid 
line). The location of the clusters are marked with a star '★' sym- 
bol. Bottom window: differences between the relative spatial 
densities using the correct (Eq. |2) and approximate (Eq. [3]) ex- 
pressions. 



Region 


d 
(pc) 


/ 

(deg) 


b 

(deg) 


Upper Scorpius 


145 


350 


+20 


a Orionis 


385 


207 


-17 


Pleiades 


130 


167 


-24 




1 2 3 

l-Ks [mag] 



spatial densities of stars and brown dwarfs of spectral type i in 
the thin disc, n{d)/no, using the correct (Eq.|2]) and approximate 
(Eq. [3) expressions are not larger than 0.2 % in the direction 
of three representative open clusters, whose galactic coordinates 
and distances are given in Table[2] The differences are two orders 
of magnitude lower (< 0.002 %) at heliocentric distances d < 
lOOpc, where the expected number of late-type dwarf contami- 
nants is maximum. These differences are very small in compari- 
son with the errors in the determination of the main parameters 
of the Galaxy disc model and, as we will see later, the spatial 
densities for each spectral type. Besides, the linear R(d, I, b) ap- 
proximation enormously simplifies next analysis. 

All n(c/)/«o curves in the top panel of Fig. [TJ show the same 
local density at null heliocentric distance, which is for an object 
of spectral type i: 

ni(R = R e ,Z = Z Q ) = n lu e = n eJ . (8) 

Upper Scorpius is in the Galactic north hemisphere, like the Sun, 
so the line of sight of the search does not cross the Galactic 
plane. Thus, in contrast to the Pleiades and cr Orionis, the Upper 
Scorpius relative star density has not a maximum at Z = (and 
nAj{d, b) and dsjid, I, b) do not change sign). For the Pleiades 
and <x Orionis, the peak of n(d)/na at dsinb = -Z Q (when the 
line of sight intersects the Galactic plane) is not exactly 1.00, 
but there is an error of about 3 % because the galactic coordi- 
nate system "is necessarily based on a galactic plane through 
the Sun", which is the H i principal plane (Blaauw et al. 1960). 
Finally, the spatial density in the direction to Upper Scorpius is 
smaller than in the direction to the other clusters at relatively 
small heliocentric distances, but it gets larger at more than about 
250 pc because the OB association is roughly in the direction of 
the Galactic Centre. 



Fig. 2. Simulated /ssa vs. /ssa - ^s,2mass colour-magnitude dia- 
gram of a hypothetical cluster survey. (Red) filled stars (10): hy- 
pothetical confirmed cluster members with spectroscopic signa- 
tures of youth; (blue) open stars (3): hypothetical cluster mem- 
ber candidates based on photometry; (blue) small dots: probable 
fore- and background sources; (red) thick dot-dashed line: se- 
lection criterion used for separating cluster members and candi- 
dates from probable fore- and background sources; (black) thin 
solid/dotted lines: "strips" where the number of contaminants 
are computed. The height of the strips were set to 1 .5 mag for 
clarity. Small dots and open stars are actual data from a multi- 
band (BRIssaJHKs,2mass) 90 arcmin-radius survey centred on 
HD 221356 (aka K6 3 A), an F8V star with an M-L-type binary 
companion in a very wide orbit (Gizis et al. 2000; Caballero 
2007). The mean sequence of the hypothetical cluster is paral- 
lel to the thick dot-dashed line, shifted 0.5 mag to the red. Filled 
stars were randomly generated by adding a 0.5 mag-amplitude 
noise to the mean cluster sequence (i.e., the selection crite- 
rion is ~ 1 cr shifted to the blue with respect to the mean clus- 
ter sequence). [SSA: SuperCOSMOS Science Archive (Hambly 
et al. 2001); 2MASS: Two-Micron All Sky Survey (Skrutskie 
etal. 2006)]. 



2.3. Integrated number of contaminants 
2.3.1. The pseudo-algorithm 

As introduced in Section QJ the number of possible field-dwarf 
contaminants in a survey can be computed from the correspond- 
ing colour-magnitude diagram. In Fig. [2] we illustrate the com- 
putation of contaminants based on a simulated colour-magnitude 
diagram. The pseudo-algorithm of integration is identical for any 
combination of "red" (e.g., K s ) and "blue" (e.g., /) passbands. 
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The general procedure to compute the integrated number of 
field-dwarf interlopers consists on: 

1. Defining horizontal "strips" in the colour-magnitude dia- 
gram. They are discrete rectangular regions in the diagram 
that broke the area of selected sources in pieces for a suit- 
able integration. The blue border of each strip is defined by 
the selection criterion; the choice of the red border is, how- 
ever, rather subjective, because observers tend to select as 
cluster member candidates all the sources to the red of the 
selection criterion. A reasonable choice of the red border is 
the colour of the reddest source in each magnitude interval, 
but it can be extended up to very red colours (indicated with 
dotted lines in Fig. [2j». The height of the strips (i.e., the size 
of the integration step) depends on the required accuracy; a 
width of ~0.5 mag can be a compromise between a fast com- 
putation and fidelity in the reproduction of the shape of the 
selection criterion (in this case, a straight line). This height 
cannot be less than the typical photometric uncertainty. 

2. Determining the spectral types of the possible interlopers in 
each strip. They are derived from the colours at the strip bor- 
ders and from a spectral type-colour relation. See below for 
caveats on this determination. 

3a. Setting the minimum and maximum heliocentric distances d\ 
and di at which an interloper of a given spectral type would 
lie within each strip, based on the corresponding absolute 
magnitudes, M,\. The computation of M,\ must account for 
the apparent magnitude m^, the heliocentric distance d and 
the interstellar extinction A,\ (the error in the Mx -spectral 
type relation is considered in item #2). 

3b. Setting the heliocentric distance d* at which the line of sight 
crosses the Galactic plane in case the minimum and max- 
imum heliocentric distances are in different Galactic hemi- 
spheres. 

4. Integrating the spatial density between the minimum and 
maximum heliocentric distances to get the number of inter- 
lopers of a given spectral type per surface unit in each strip. 

5. Summing the number of interlopers of all the considered 
spectral types and in all the strips. 

6. Multiplying the value by the survey area, S, to get the inte- 
grated number of contaminants. 

Regarding item #2, a field dwarf with intrinsic colour bluer 
than the selection limit can still lie in the selected region of 
the colour-magnitude diagram because of photometric error. 
Moreover, a spectral type does not correspond to one absolute 
magnitude and one colour value but to a range of magnitudes 
and colours. It is necessary to take these effects into account 
when calculating the contamination, otherwise the result will 
be underestimated. One way for accounting for the photometric 
scattering of the data (including the broadening of the main se- 
quence and the survey photometric uncertainties) may be choos- 
ing one earlier spectral subtype of the bluer interloper than ex- 
pected from the selection criterion. At the late considered spec- 
tral types, a variation of one subtype corresponds to an average 
colour variation of, e.g., A(7 - J) at 0.2-0.3 mag (see below), 
which may be a too large photometric uncertainty for a typi- 
cal survey. An alternative solution is simply to shift bluewards 
the selection criterion actually used in the colour-magnitude di- 
agram by certain amount (~0.05-0. lOmag). This amount would 
depend on the survey, magnitude range, and considered spectral 
types. 

Once the main parameters (TableQ]) and the exponential form 
of the Galaxy thin-disc model (Eq. [T]) are known, as well as the 
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spatial densities, colours and absolute magnitudes of the possible 
contaminants as a function of spectral type (Section 12. 4K the 
integration can be carried out. 

2.3.2. The integration 

In this Section we provide the necessary tools to integrate the ex- 
ponential spatial density in the linear approximation (Eq.|5j. In a 
typical photometric search program, assuming that such a search 
does not span an excessively wide area (i.e., smaller than a few 
tens of square degrees), the search volume can be approximated 
as a cone with an apex length equal to the heliocentric distance 
of a given source detectable within the photometric depth, and a 
base corresponding to the projected area on the sky of the search, 
S. The axis of the cone is defined by the line of sight (i.e., the 
galactic latitude and longitude of the centre of survey area). The 
use of a cone makes the computation of late-type contaminants 
much easier, since it is not necessary to know the actual shape of 
the survey. 

The number of late-type dwarfs of spectral type i that con- 
taminate the survey comes from the integration in a truncated 
cone of volume V with a non-uniform density n(r): 

Nt = / m(d)dV = I JT I m(z, <P, a)dV, (9) 

where the volume differential in the coordinate system of the 
cone is z ~*™" dzdipda (being z the regular z-coordinate, <p the 
angle of a projected vector with the positive X-axis, measured 
counterclockwise, and a the angle between the vertical axis, the 
cone vertex and any other point of the cone). When the galactic 
latitude and longitude of the field of view are known, the density 
n(r) depends only on the heliocentric distance, d (the "height" 
of the cone). Using the simplified expression of the density of 
thin-disc dwarfs (Eq.|5j, the integration between the heliocentric 
distances d\ and d 2 provides: 

» t rdi , . t j r ln a j. rarctan (fe/A) sinQ , 

Ni = J di ni (z)z 2 dzj Q d0j Q z^da = 

= tt(!) 2 J^n^-^dz, (10) 

In the pyramid-cone approximation of height h and base radius 

b, the value ^(|) is simply the angular area of the search, S, 

which must be expressed in rad 2 . If d\ and d2 are each in a 
Galactic hemisphere (i.e., the line of sight crosses the Galactic 
plane), then the integral is detached into two parts: 

N t = s(n A+l jf e^z 2 dz + n A __ ^ e~^ z 2 dz) , (11) 

where d„ is the heliocentric distance at which n(d) has a max- 
imum (Z = 0) and the n A and d B auxiliar variables change be- 
haviour. The exponential integral has a primitive, which is: 

j£ e -^z 2 dz = [(-l)e'i(d B z 2 + 2d 2 B z + 2d 3 B )] d J = (12) 

= I(d B ;di,dj). (13) 

Hence, the final expression for the total number of dwarfs in a 
truncated cone with the minor and major bases at heliocentric 
distances d\ and da, respectively, is: 

Niin, r 2 ) = S [n A+J I(d B+ ; d\,d t ) + n A _ t I(d B J,d», d 2 )\ . (14) 
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M5 LO L5 TO T5 



Sp. type 

Fig. 3. Local spatial density of ultracool dwarfs as a function of 
spectral type from the data in Table [3] 

2.4. Possible late-type contaminants 

We have studied the contribution to contamination in photomet- 
ric searches for very red sources by field dwarfs with spectral 
types from M3 V to T8 V. This election is based, on the one hand, 
on the earliest spectral type of a brown dwarf in a ~1 Ma-old 
star-forming region, which is about M5, and the large-amplitude 
photometric variability that can be present in the brightest young 
substellar objects, which can be as large as 0.7 mag (Caballero 
et al. 2006). The maximum colour variation observed could 
be even larger than the amplitude of photometric variability. 
Because of that, if the surveys in each passband are not simulta- 
neous, very young variable brown dwarfs can display colours of 
objects a few sub-types earlier than its actual spectral type. On 
the other hand, the late-type bound is set by the latest-type brown 
dwarfs observed to date (e.g., Burgasser et al. 2000; Tinney 
et al. 2005; Looper, Kirkpatrick & Burgasser 2007; Lodieu et al. 
2007a; Warren et al. 2007). The total number of contaminants 
in the survey will be the sum of the number of field dwarfs of 
spectral type i in the truncated cones of search: 

Num. 

N=^ J N i (r u ,r x d, (15) 
;=i 

(i =M3^1V, M4-5V, M5-6V... T7-8V). 

The Johnson /-band absolute magnitudes, I - J and J - K s 
colours (J and K s in the 2MASS system) and local spatial densi- 
ties from the literature of the 25 considered intervals of spectral 
types are given in Table [3] For the M3-L1 dwarfs, we used the 
tabulated values from West et al. (2008) of /sdss - J, then cal- 
culated the / - /sdss corrections using the template spectra of 
Bochanski et al. (2007) to get the I - J colours. Absolute /-band 
magnitudes were computed from absolute /-band magnitudes 
compiled by Hawley et al. (2002) and using our I - J values. 
J - K s colours are also from West et al. (2008). For the L8-T0 
dwarfs, J-K s values are from Dwarf Archives . org, taking av- 
erage of objects with uncertainties 5(J - K s ) < 0. 1 mag and not 
known to be binary. The absolute / and I - J colours are based 
on the polynomial relations of Tinney, Burgasser & Kirkpatrick 
(2003 - note the backtrack ml- J colours). The remaining IJK S 
photometry has been taken from Dahn et al. (2002) and Vrba 
et al. (2004). 



Table 3. Characteristics of possible late-type field dwarf con- 
taminants from the literature . 



Spectral M, I - J J - A* s n Q d b I=25 

type (mag) (mag) (mag) (1(T 3 pc~ 3 ) (pc) 



M3-4 V 


9.07 


1.83 


0.84 


19.6 


7700 


M4-5 V 


10.14 


1.80 


0.86 


16.3 


5700 


Mj- 6 V 


1 i.y 1 


i /n 
2.4 / 


n on 
0.90 


n o 1 

y.si 


1 1 nn 


AJt£. H \ 7 

Mo—/ V 


1 1 OA 

I2.yu 


1 71 

2. 11 


n nQ 
0.93 


a nn 

4.yu 


linn 
11UV 


A/fH O \ 7 

Ml— 8 V 


14. to 


i i/i 
3.24 


n no 
0.98 


1 nm 

i.yu/ 


1 Qnn 
1 JUU 


A to n \ 7 

M8-9 v 


14.68 


3.54 


1 n£ 

1.06 


1.285 


1 100 


My— LU V 


1 3.23 


i on 
3.80 


1 11 
1.12 


0. /68 


o/i n 
841) 


t n 1 ^ 7 
LU— 1 V 


1 C n 

1 J. 11 


3.55 


1 1"7 

1.2/ 


n t^m 
U.03 / 


oin 
821) 


T 1 1 \7 

Li— 2 V 


1 C £Ci 

I J. W 


i cn 
3.50 


1 Q O 

1.38 


U.861 


/on 
08L) 


T 1 1 \ 7 

LI— 5 V 


1 H 1 Ct 

lo. iy 


3.02 


1 Al 

1.43 


1 nnn 
l.UUU 


550 


T 1 A \ 7 

Li— 4 V 


16.6/ 


3. 11 


1 An 

i.4y 


n OQ A 

0.834 


43L) 


L4-5 V 


17.14 


3.81 


1.54 


0.819 


360 


L5-6V 


17.61 


3.90 


1.60 


0.869 


290 


L6-7V 


18.07 


4.06 


1.66 


0.785 


240 


L7-8V 


18.52 


4.12 


1.71 


0.644 


190 


L8-9V 


18.90 


4.22 


1.81 


0.462 


160 


L9-T0V 


18.95 


4.40 


1.73 


0.308 


160 


T0-1 V 


19.01 


4.63 


1.58 


0.220 


160 


T1-2V 


19.04 


4.94 


1.17 


0.241 


150 


T2-3 V 


19.04 


5.19 


0.81 


0.406 


150 


T3^1V 


19.07 


5.23 


0.51 


0.788 


150 


T4-5 V 


19.22 


5.01 


0.29 


1.42 


140 


T5-6V 


19.61 


4.76 


0.11 


2.38 


120 


T6-7V 


20.35 


5.25 


0.00 


3.45 


85 


T7-8V 


21.60 


5.50 


-0.06 


4.82 


48 



" The sources of the M ; absolute magnitude, I- J and J - K s colours, 
and spatial densities are obtained or derived by us from data in 
the literature (Kirkpatrick et al. 1994; Dahn et al. 2002; Hawley 
et al. 2002; Cruz et al. 2003; Tinney et al. 2003; Vrba et al. 2004; 
Bochanski et al. 2007; Burgasser 2007; West et al. 2008). See the 
text for further details. 

b di = 25 is the heliocentric distance of a dwarf with an apparent magni- 
tude / = 25.0 mag, assuming an interstellar extinction proportional 
to the heliocentric distance. At the end of Section [2T4l there is a brief 
discussion on this assumption. 



The local spatial densities (n G ) and the spatial densities at 
R = R and Z = («o) are related through Eq. [8] A pictori- 
cal representation of n versus spectral type is given in Fig. [3] 
Spatial densities of M3-7 and M7-L0 dwarfs are measured val- 
ues from Kirkpatrick et al. (1994) and Cruz et al. (2003), re- 
spectively. For L0-T8 dwarfs, we used predicted values from 
Burgasser (2007) based on a Monte Carlo simulation assuming 
a mass function dN/dM oc M~ a5 (similar to estimates from var- 
ious star forming regions - e.g., Luhman et al. 2000; Caballero 
et al. 2007), a normalisation of 0.0037 pc~ 3 over the range 0.09- 
O.1OM (Reidetal. 1997), a binary fraction of 10% and a binary 
mass ratio distribution of P{q) oc q 4 (where q = M2/M1). These 
values are consistent with empirical space density estimates by 
Burgasser et al. (2002), Cruz et al. (2007) and Metchev et al. 
(2008), and predicted densities by Deacon & Hambly (2006), 
although they differ significantly from predictions by Nakajima 
(2005) for L5-T7 dwarfs. 

Using the expression I - Mi = 5 log d — 5 + Ai, a deep sur- 
vey with a limiting magnitude I\\ m = 25.0 mag would be able to 
detect a T5-6 dwarf at an heliocentric distance d ~ 120pc and 
an M5-6 dwarf at about 3.1 kpc to the Sun (see last column in 
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Table [3}. In the direction of the Pleiades, such an hypothetical 
M-type "field" dwarf would be located at about 1.9kpc below 
the Galactic plane at almost z = 6 x h z . In this location, the 
spatial density of thin-disc objects of any spectral type is about 
10~ 3 times the local density n Q - It is evident from this example 
the necessity of the integration along the cone of search using a 
non-uniform density. Following this assumption, we should also 
use an interstellar extinction that does not linearly increase with 
heliocentric distance. Instead, it should be proportional to the ac- 
tual interstellar dust content. In our model, we assume that the 
Galactic dust content has much wider length and height scales 
than for stars. 

3. A practical application: L and T dwarf 
contamination in a Orionis 

3.1. Preliminaries 

The cr Orionis cluster (~3Ma, ~385pc) contains one of the 
best known substellar populations, from the hydrogen burning 
mass limit down to a few Jupiter masses (see a bibliographic re- 
view in Caballero 2008b). It could be the most favourable site 
known to date for studying the opacity mass limit for forma- 
tion of objects via fragmentation in molecular clouds, that may 
lie just below 0.003 M (Rees 1976; Bate, Bonnell & Bromm 
2003 and references therein). Even accounting for the L-type 
cluster members and candidates in the Pleiades (Martin et al. 
1998; Dobbie et al. 2002b; Moraux et al. 2003; Bihain et al. 
2006), Chamaeleon and Lupus (Allers et al. 2006; Jayawardhana 
& Ivanov 2006), <x Orionis contains by far the largest number of 
L-type objects with spectroscopy (12 - Zapatero Osorio et al. 
1999, 2000; Martin et al. 2001; Barrado y Navascues et al. 2001, 
2003). While there is reasonable criticism on the cluster mem- 
bership of some of the faintest objects (e.g., S Ori 47, LI. 5+0. 5 
- McGovern et al. 2004), severeal of them are known to be ex- 
tremely young based on near-infrared excess and extraordinary 
Ha emission of up to -700 A (S Ori 55, S Ori 71 - Barrado y 
Navascues et al. 2002; Zapatero Osorio et al. 2002a; Caballero 
et al. 2007) or flux excess in the IRAC Spitzer bands (Zapatero 
Osorio et al. 2007b; Scholz & Jayawardhana 2007). Besides, 
cr Orionis also possesses the only extremely young T-type cluster 
member candidate identified so far (S Ori 70; it was introduced 
in Section [TJ- Casewell et al. 2007 and Bouvier et al. 2008 have 
recently identified T dwarf candidates and cluster members in 
the Pleiades [-120 Ma] and the Hyades [-600 Ma]). 

The contamination by field late-M dwarfs in the cluster is 
well bounded. Caballero et al. (2007) identified 30 objects as 
bona fide <x Orionis members with features of extreme youth 
amongst 49 candidate cluster members fainter than I = 16 mag 
selected from an I vs. I - J diagram. They also reported the 
detection of two L-type field dwarf candidates in the direction 
of the cluster, even fainter (and much more distant) than Cuby 
et al. (1999)'s T dwarf. The ratio between very low-mass con- 
firmed and candidate cluster members in cr Orionis is compa- 
rable to that in the Orion Nebula Cluster (Lucas et al. 2001; 
Slesnick, Hillenbrand & Carpenter, 2004), but much larger than 
in other older open clusters like the Praesepe (Pinfield et al. 
1997, 2003; Adams et al. 2002; Gonzalez-Garcia et al. 2006) or 
the Hyades (Reid & Hawley 1999; Gizis, Reid & Monet 1999; 
Dobbie et al. 2002a). In any case, some of the faintest substellar 
member candidates in cr Orionis have no astrometric confirma- 
tion or irrefutable youth signatures. Therefore, while there are 
no new data on cluster membership, it is necessary to accurately 
determine the number of possible fore- and background contami- 



nants among selected photometric cluster member candidates for 
investigating the cluster Initial Mass Function. We have applied 
our Galaxy and late-type dwarf data and mathematical tools to 
estimate the number of fore- and background L and T dwarfs as 
a function of depth in a hypothetical survey area of 1 deg 2 (S « 
0.000 30rad 2 ) in the cr Orionis cluster. 

We have assumed that the cluster member selection is ex- 
clusively based on I - X colours (Johnson 7), being X a near- 
infrared passband (e.g., JHK S ). The X-band observations are 
deep enough to detect within the completeness all the L- or T- 
type object identified in the /-band images^. The combination 
of a red optical filter and a near-infrared one is an effective way 
(i.e., in telescope time, in simplicity) for selecting very faint red 
objects in young open clusters. Although an / - /T s -based search 
would be better for identifying late-L-type objects, the increase 
in thermal background at 2 /jm of near-infrared detectors favours 
the use of the J passband, especially for T-type objects (that have 
quite blue J - K s colours); a compromise might be found at an 
I - H survey. Likewise, imaging at passbands bluewards of / 
is prohibitive due to the extremely dimness of ultracool cluster 
member candidates at such wavelenghts (e.g., a cluster member 
with / = 24 mag is expected to have R « 27-28 mag). Finally, the 
number of interlopers with pure optical (R - 1, 1 — Z) and pure 
infrared {J - K s ) passbands is larger than with a combination of 
red optical and near-infrared filters, because of the larger slope 
of the spectral energy distribution of young late-type objects at 
the optical/infrared boundary (i.e., redder I — J, I - K s colours). 

As explained in Section l2.3.1l one of the steps of the integra- 
tion of the number of contaminants is setting the minimum and 
maximum heliocentric distances d\ and dz specified in Eq. [14] 
For the computation of the heliocentric distance of an object of 
apparent and absolute magnitudes / and M/, one must find the 
root of the smooth function f(d) = 5 log d — 5 + Aj — (/ - M/), 
where the extinction A/ = 0.482Ay is proportional to the dis- 
tance (A/ = ad) and A v = 3.09E(B - V) (Rieke & Lebofski 
1985). The constant a depends on the colour excess E(B - V)* 
and distance d+ towards cr Orionis (e.g., Lee 1968; Caballero 
2008a). Quantitatively: E(B - V) = E(B - V)*(d/d*). To find the 
root, we have used the Newton method {d n+ \ = d„ - f(ri)/ f'(ri), 
n = 0, 1, 2...), that quickly converges in 3-6 iterations with a 
suitable initial value of do- 

In Table |4] we summarize our results and provide the num- 
ber of L0-T8V-type dwarfs in 1.0 mag-width "strips". As dis- 
cussed below, not all the late-type field dwarfs in the survey 
area are true contaminants (e.g., L-type dwarfs do not contam- 
inate the cr Orionis / vs. / - J colour-magnitude diagram at 
/ > 24.5 mag). The uncertainties come from the errors in the 
parameters of the thin disc (Table[TJ, the cluster heliocentric dis- 
tance (d = 330—440 pc), the spectral type-absolute magnitude 
relation, the colour excess, and the local spatial densities. We 
use conservative errors of 0.2 mag and 0.02 mag for the spec- 
tral type-M/ relation and the interstellar extinction, respectively. 
Variations in the solar galactocentric distance R Q barely affect 
the results. The most important contributors of error are the ra- 
dial scale length h r , that is the Galactic parameter with the largest 
intrinsic uncertainty (> 40 %), and the local spatial densities. The 
mass function index that Burgasser (2007) used for deriving n 
(a = +0.5; dN/dM oc M~ a ) may be too large according to re- 
cent observational results (Lodieu et al. 2007b; Metchev et al. 



1 For example, the completeness magnitudes in the representative 
survey by Caballero et al. (2007) were / comp i. = 23.3 mag, / comp i. = 
20.6 mag, allowing the detection of all the objects within the complete- 
ness in the optical redder than I - J = 2.7 mag (~M7V). 
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Table 4. Number of expected field L- and T-type dwarfs in 
a 1 deg 2 -area survey towards <x Orionis as a function of the 
Johnson /-band magnitude and spectral type interval". 



Al 


L0-5V 


L5-T0V 


T0-5V T5-8V 


(mag) 









21.0- 


-22.0 


i r - 5 

-0.4 
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" Poissonian localized overdensities are not taken into account. 



2008). A flatter substellar mass function (a * 0.0) would lead 
to lower local spatial densities and, therefore, lower contami- 
nation rates. However, the choice for Burgasser (2007)'s single 
power-law mass function (from where L0-T8V local densities 
in Table Q] were obtained) seems more reasonable to us, since it 
explains the raise of the mass function found by many other au- 
thors in the planetary regime (see Section Q}. Besides, Metchev 
et al. (2008)'s values are consistent with Burgasser (2007)'s ones 
within the 95 % confidence limits. Maximum differences with 
other determinations of substellar densities in the literature are 
of the order of 20 %. Our estimation of the uncertainties in n Q 
is, however, more conservative: maximum and minimum sur- 
face densities come from using very different indices of the lo- 
cal mass function, a — and +1 (see table 5 in Burgasser 2007). 
The error in n strongly depends on the spectral type, and ranges 
from ~ 12 % at L0-1 V to -50 % at T7-8 V. 

Our computations surpass and complement the estimate of 
contamination in the cr Orionis by Caballero et al. (2007). The 
number of interlopers shown in Table|4]are consistent with those 
given in Caballero et al. (2007) accounting for the different sur- 
vey area, input absolute magnitudes M/ and local densities n Q 
(especially in the L5-T7V interval), and width of strips (i.e., 
I- J colour). They estimated that ~4 L-type dwarf contaminants 
(~2 early, ~2 late) populated their least-massive bin of the Initial 
Mass Function, that contained 1 1 planetary-mass object candi- 
dates (i.e., the contamination rate is 36 %); our results support 
those calculations. A coarse extrapolation towards fainter mag- 
nitudes of the surface densities provided by Burgasser (2007) 
also matches our results (he did not account, however, for the 
Galactic structure). Similar values as in Table [4] are expected 
for 1 deg 2 -wide surveys of any kind (young open clusters, extra- 
galactic) at intermediate galactic and high latitudes of the same 
magnitude depth. 

3.2. Discussion 

In short, from Table |4] the number of ultracool dwarfs in deep 
surveys is larger for earlier spectral types and for the faintest 
magnitude intervals. The contamination by T dwarfs is very low 
for survey depths / < 24 mag. Surveys deeper than this value 
are strongly affected by L- and T-type dwarf interlopers. From 
Fig. [4] it is not clear what happens to the cr Orionis luminos- 
ity function at / > 24 mag. Currently, the deepest optical sur- 



veys in young clusters in general, and in cr Orionis in particu- 
lar, reach limiting magnitudes / = 24-25 mag (Caballero et al. 
2007 and references therein - see also Comeron & Claes 2004 
for a very deep near-infrared survey in Chamaeleon I), although 
some ultra-deep extragalactic surveys can go 2-Anmg fainter. 
The magnitudes (/ * 22.8 and 23.2 mag) and measured spec- 
tral types (L5. 0+1.0 and L3. 5+2.0 - Barrado y Navascues et al. 
2001; Martin et a al. 2001) of the faintest cluster objects with 
unambiguous youth features, S Ori 65 and S Ori 60 (they display 
flux excess at 8.0 and/or 5.8 fim - Zapatero Osorio et al. 2007b; 
Scholz & Jayawardhana 2008) suggest that the L-T transition in 
the cluster occurs at about / ~ 23.5 mag (field) T0-1V dwarfs 
are about 1 .5 mag fainter than L5-6 V ones - however, this as- 
sumption might not be necessarily right for low-gravity sources 
due to pressure effects on alkali line absorption). S Ori 69, the 
second faintest cluster member candidate {I a 23.9 mag) has a 
tentative TO: spectral type from a near-infrared spectrum (Martin 
et al. 2001), which supports our estimate (the semi-colon indi- 
cates uncertainty in the classification). Therefore, fainter clus- 
ter members are expected to have T spectral types, with colours 
I-J> 4.5 mag and most probable masses derived from state-of- 
the-art theoretical models at M < 0.006 M (we do not discuss 
in this paper the cluster membership status of S Ori 70). 

The actual number of cr Orionis cluster members and can- 
didates of L0-5 and L5-T0 spectral types at the correspond- 
ing magnitudes / ~ 21.0-22.5 and 22.5-24.0mag in a sur- 
vey is larger than the expected number of dwarf contaminants 
of the same spectral type (e.g., Zapatero Osorio et al. 2000). 
The difference, although significant, requires a careful spectro- 
scopic/photometric/astrometric follow-up because fluctuations 
about the mean could mask the data. A different space den- 
sity of contaminants towards <x Orionis by some Poissonian fac- 
tor would easily lead the L-type contamination rate to grow up 
from 36 to ~70% of more. Localized overdensities of dwarfs 
in the Galactic thin disc are not unusual (e.g., Juric et al. 2008 
and references therein). Field L dwarfs stop contaminating the 
cr Orionis colour-magnitude diagram at / w 24.5 mag (i.e., 
I - J > 4.5 mag). Since the interstellar reddening is apprecia- 
ble only at very large heliocentric distances (A/ = ad, being 
a « 2 10~ 4 mag pc~'), where the spatial density is extremely low, 
contamination by M dwarfs is not expected, either. 

Surveys with completeness magnitudes 7 comp i. « 27-28 mag 
should be required to search for the opacity mass limit (at M < 
0.003 M Q ), and field T dwarfs are expected to be the major 
colour-selected contaminants at these magnitudes. While the 
contamination by T dwarfs in the brightest magnitude inter- 
vals in Table [4] is very low with respect to the contamination 
by L dwarfs (by a factor ~100), it dramatically increases in 
the faintest magnitude intervals. For example, for Al = 28.0- 
29.0 mag, the number of L0-L5V dwarf contaminants reaches 
the maximum limit of ~300 objects (i.e., the survey is so deep 
that is able to detect almost all the early L dwarfs of the Galaxy 
in that direction), whereas there are only 3.6 times less T dwarfs 
than L dwarfs for the same Al. From the data in Table [4] we ex- 
pect about ~45 T-type dwarfs in the foreground of cr Orionis 
in an ultra-deep (/ w 27-28 mag) survey. However, the most 
optimistic estimations of a large index a = +0.6 of the cluster 
mass spectrum (Caballero et al. 2007) extrapolated towards only 
0.001 M Q predict <30 planetary-mass objects with T spectral 
types in cr Orionis (Zapatero Osorio, priv. comm.). Assuming 
less favourable conditions, we forecast no more than ~10 T- 
type cluster members. These values depend on assumptions that 
require confirmation: e.g., the distance and age of the cluster, 
the actual value of the slope of the mass spectrum at Jovian 
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Fig. 4. I vs. I - J colour-magnitude diagram of the cr Orionis 
cluster. Big filled circles with error bars: cluster members and 
candidates (the red outlier is the T-type object S Ori 70); solid 
and dashed lines: COND03 and DUSTY00 3 Ma-old isochrones, 
respectively, shifted at d — 385 pc (Chabrier et al. 2000; Baraffe 
et al. 2003); dotted lines: dwarf sequence in Table [3] shifted at 
distance moduli m — M= 3.0, 6.0, 9.0, 12.0, and 15.0 mag, from 
top to bottom. The faintest data point of the COND03 theoret- 
ical isochrone (7 - / w 3.0 mag, 7 =s 29.0 mag) corresponds to 
an l.OMjup-mass object. L and T field dwarfs have approximate 
colours redder than I - J « 3.5 and 4.5 mag. 



masses, or the validity of theoretical evolutionary tracks of low- 
mass stars and brown dwarfs at very young ages and low masses 
(Baraffe et al. 2002). For example, Fig. g] shows that COND03 
models (Baraffe et al. 2003), which are especially useful for the 
T-type domain in the field, predict rather blue I - J colours for 
the least massive hypothetical cr Orionis members (with masses 
down to ~1.0Mj up ). Accurate theoretical modeling of low grav- 
ity, low temperature spectra is still work in progress. 

Even with the largest current ground facilities (2 x 10.0 m 
Keck Observatory, 10.4 m Gran Telescopio Canarias), the neces- 
sary follow-up of an ultra-deep survey for studying the opacity 
mass limit in cr Orionis would require a tremendous, prohibitive, 
observational effort. For example, the low-resolution near- 
infrared spectrum of the ~T6 object S Ori 70 (/ = 25.0±0.4mag, 
J = 19.98+0.06 mag) needed a total exposure time of 4 800 s 
with NIRSPEC at the Keck II (Martin & Zapatero Osorio 2003 
- it had, besides, a low signal-to-noise ratio). Sources with 7 ss 
27-28 mag would need correspondingly larger exposure times, 
of up to a whole night. Therefore, spectroscopy of roughly half 
of the selected cluster member candidates in such an ultra-deep 
survey (> 20 targets) would demand about one month of observ- 
ing time at the Keck II. Shallower surveys (7 < 26 mag) and 
their corresponding follow-ups are, however, much less affected 
by T-dwarf contamination (~5 interlopers in a 1 deg 2 -wide sur- 
vey) and would demand a quite more reasonable amount of time 
to be accomplished. The photometric selection becomes more 
stringent and may lower the contamination level when combin- 
ing several colour-magnitude diagrams as well as colour-colour 
diagrams and CFLi-on and Oj(f photometry. A survey with at least 



three broad filters (e.g., 7/7f s or z'77[4.5]2 an d the two narrow 
methane filters, although it multiplies the total exposure time by 
a factor ~5-6 with respect to an 77 survey, may represent an 
intermediate solution. 

4. Summary 

Late-type field dwarfs are the most important contributors to 
contamination among faint, young, cluster member candidates 
selected from colour-magnitude diagrams. The number of such 
interlopers in deep photometric surveys in clusters is of the max- 
imum importance for the study of the substellar Initial Mass 
Function. We provide a pseudo-algorithm, an expression for the 
integrated number of field dwarfs at certain heliocentric distance 
and galactic coordinates (assuming an accurate linear approxi- 
mation in the structure of the Galactic thin disc), and the abso- 
lute magnitudes, colours and local spatial densities of M3-T8 
dwarfs to compute the number of field late-type contaminants in 
deep surveys excluding the Galactic plane. 

We have applied our tools and data to an hypothetical ultra- 
deep (7 comp i = 29 mag) survey towards the young (~3Ma) 
cr Orionis cluster. We predict a rather low contamination rate 
of L and T field dwarfs with an appropriately defined selection 
criterion up to 7 = 25-26 mag. The number of contaminants at 
fainter magnitudes, where the opacity mass limit is expected to 
lie, is however very large and may preclude photometric, spec- 
troscopic and astrometric follow-up in reasonable amounts of 
time with current facilities. 

The enhancement of the contamination by T-type dwarfs in 
ultra- and very deep surveys is not only important for <x Orionis, 
in particular, and star-forming regions, in general, but also for 
extragalactic surveys of z > 6 quasars, whose colours resemble 
those of T dwarfs. 
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